Interacting particles at a metal-insulator transition 
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We study the influence of many-particle interaction in a system which, in the single particle case, 
exhibits a metal-insulator transition induced by a finite amount of onsite pontential fluctuations. 
Thereby, we consider the problem of interacting particles in the one-dimensional quasiperiodic 
Aubry-Andre chain. We employ the density-matrix renormalization scheme to investigate the fi- 
nite particle density situation. In the case of incommensurate densities, the expected transition 
from the single-particle analysis is reproduced. Generally speaking, interaction does not alter the 
incommensurate transition. For commensurate densities, we map out the entire phase diagram and 
find that the transition into a metallic state occurs for attractive interactions and infinite small 
fluctuations — in contrast to the case of incommensurate densities. Our results for commensurate 
densities also show agreement with a recent analytic renormalization group approach. 
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I. INTRODUCTION 

The metal-insulator transition (MIT) in disordered 
electronic systems has been the subject of intense re- 
search activities over the last two decades and still con- 
tinues to attract much attention. For non-interacting 
electrons in disordered systems Q] the scaling hypoth- 
esis of localization ||] can successfully predict many of 
the universal features of the MIT. However, the influence 
of many-particle interactions on the MIT is not equally 
well understood ^ and recent investigations of an ap- 
parent MIT in two-dimensional systems even question 
the main assumptions of the scaling hypothesis . In 
the single-particle case, the one-dimensional quasiperi- 
odic Aubry-Andre model is known rigorously to exhibit 
an MIT for all states in the spectrum as a function of 
the quasiperiodic potential strength fi [Q. The ground 
state wave function is extended for /i < 1 and local- 
ized for /i > 1. The system at /Uc = 1 is critical: there 
the wave functions decrease algebraically, not exponen- 
tially as in the localized case. This behavior is contrary 
to the localization of non-interacting particles in a one- 
dimensional random potential, where the ground state 
wave function becomes localized for infinite small disor- 
der strength. An ingeneous theoretical approach to the 
interplay of interactions and disorder is based on the two- 
interacting-particles (TIP) problem in one-dimensional 
random potential p|-[lT[|. Furthermore, numerical results 
for spinless fermions in a random potential at finite par- 
ticle density have given additional insight [|2 13|. In gen- 
eral, these investigations have shown that changes in the 
wave function interferences due to many-particle inter- 
actions can lead to a rather large enhancement 
of the localization lengths in one and two dimensions 
| p^ , p^p7| . Recently, we examined the TIP problem in a 
quasiperiodic potential by means of the transfer-matrix 
and the decimation method together with a careful finite- 



size-scaling analysis [^^Q . We found that it exhibits the 
MIT at /ic = 1 as in the single-particle case — indepen- 
dent of interaction. 

As an independent extension of these low-density TIP 
results, Chaves and Satija have studied a model of 
nearest-neighbor interacting spinless fermions ||2l|-p^ at 
finite particle density in the same quasiperiodic potential 
by means of Lanczos diagonalization for small systems 
up to chain length M — 13. They find evidence for a 
critical region in which the behavior of the charge stiff- 
ness p^-p6t is different from its behavior in the metallic 
and localized regimes. In order to reach much larger 
system sizes for interacting systems, we employ the nu- 
merical density-matrix renormalization group (DMRG) 
p7t which has been shown to be very useful . In par- 
ticular, the ground state properties of interacting many- 
particle systems in one dimension can be obtained very 
accurately ||2^,^. In the present paper, we shall study 
the quasiperiodic model of Ref. at various densities 
and interaction strengths by DMRG and compare our 
results to the TIP data. 

For the interacting many-body system in a quasiperi- 
odic potential we recover the transition at /^c = 1 in- 
dependent of interaction, provided we consider densities 
like p = 1/2 which are incommensurate compared to the 
wave vector of the quasiperiodic potential — an irra- 
tional multiple of tt. Thus, the low-density TIP case 
is comparable to finite but incommensurate densities. 
On the other hand, for commensurate densities, we find 
that the system can be completely localized even for 
/i ^ 1, due to a Peierls resonance between electronic 
and quasiperiodic potential degrees of freedom. Whereas 
for repulsive interactions the ground state remains local- 
ized, we find a region of extended states for attractive 
interaction due to the interplay between interaction and 
quasiperiodic potential. The behavior may be described 
by a weak-coupling renormalization group (RG) treat- 
ment pi]]. Thus, the physics of the model is dominated 
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by whether the density is commensurate or incommensu- 
rate and only then by interaction effects. 

The paper is organized as follows. The finite-density 
many-body system is introduced in section ||. Results 
for incommensurate and commensurate densities are pre- 



five lattice sweeps and keep 300 (small systems) to 500 
(larger systems, M =89 and 144) states per block. As our 
observable of the phase transition we choose the phase 
sensitivity 



sented in section III and [V , respectively. We summarize 
and conclude in section |V|. 



II. THE NEAREST-NEIGHBOR HAMILTONIAN 
AT FINITE DENSITY AND THE DMRG 

Let us consider N interacting spinless fermions on a 
ring of circumference M in the Aubry- Andre potential 
such that 

H = -tY^ (cL+iC,„ + h- c. 

ni 

n,„+in„i + 2/i^rtmCOS (am -I- /?). (1) 

m m 

The operators cj„, Cm, and rim, respectively, denote as 
usual Fermi creation, annihilation and number operators; 
2/icos(am + /3) is the quasiperiodic potential of strength 
fi with a/2TT being an irrational number. /3 is an arbi- 
trary phase shift and we choose a/2TT = {Vb — l)/2, i.e., 
the inverse of the golden mean. In addition, we set t — I. 
The particle density is p = N/M. The model has been 
studied extensively in two limits, namely, independent 
particles y = 0, and the clean case, ^ — pl|-p3|. 
The first case corresponds of course to the single-particle 
problem discussed briefly in section^. In the second case, 
the model can be mapped onto the anisotropic Heisen- 
berg (XXZ) model for which a closed Bethe-ansatz solu- 
tion exists pl|-p3{ . It shows three distinct phases at zero 
temperature. For half filling (p — 1/2) and strong repul- 
sive interaction {V > 2) the system is a charge-density- 
wave insulator. For weak and intermediate interaction 
strength and away from half filling it is a metal and can 
be described as a Luttinger liquid with linear energy dis- 
persion and gapless excitations [^,^. The metal, at 
least, is separated for all fillings by a first-order transi- 
tion at V = —2 from an insulator where the fermions 
form clusters (phase-separated state), corresponding to 
the ferromagnetic state of the spin model |2^ . 

In a previous analysis pO[ | of the interacting quasiperi- 
odic model (|l|) , a large enhancement of the Drude weight 
D (or Kohn stiffness) and the superconducting fluc- 
tuations in the ground state was found near the first- 
order transition at y « — 2 using exact Lanczos diag- 
onalization. However, the system sizes attainable by 
the diagonalization were restricted to M < 13. Using 
the DMRG, it is possible to extend the tractable sys- 
tem lengths to about M w 100 — 200. We use the finite 
lattice algorithm for non-reflectionsymmetric models as 
described in Ref. pSl. In our simulations we perform 



MAE = M(-1)^[£;(0) - E{tt)] 



(2) 



of the ground state |32l, which is connected to D in the 
clean case (Luttinger phase) and equal to n^D for non- 
interacting fermions. Here E(^) measures the reaction 
of the system due to a twist in the boundary condition, 
CM+i — exp(i$)ci |24M. The prefactor (—1)^ cancels 
the odd-even effects 113617 In addition, it is believed that 
the phase sensitivity matches the character of the wave 
function. It does not depend on system size if the wave 
function is extended, e.g., for /i = 0, and it decreases 
exponentially for large systems if the wave function is lo- 
calized [|o| . As argued in Ref. this can be transferred 
to the critical case. Therefore, we suppose that the phase 
sensitivity decreases algebraically with increasing system 
size, if the wave function is critical. Thus, we have to 
compare at the very least three different chain lengths in 
order to characterize the length dependence of the wave 
function. 



III. INCOMMENSURATE PARTICLE DENSITIES 

In order to systematically study finite-size effects, we 
should in principle compute MAE for fixed particle den- 
sity p and increasing M. However, with fixed a/2Tr = 
(V5 - l)/2, the quasi-periodic potential is incompatible 
with periodic boundary conditions and the results for 
MAE depend on the choice of /3. Averaging over many 
potentials with different /3 can in principle be used in or- 
der to obtain /3-independent results. However, we have 
found that the statistical fluctuations for densities like 
p = 1/2, 1/3, 1/4 at fixed a/27r = (\/5 - l)/2 are rather 
large such that we could not detect any significant M 
dependence of the phase sensitivity. 

As is customary in the context of quasiperiodic systems 
prf , the value of a/27r may be approximated by the ratio 
of successive Fibonacci numbers — = f„_2 + ^"n-i = 
0, 1, 2, 3, 5, 8, 13, .... In this way, choosing a/(27r) = 
Fn-i/Fn and M = F„, we can retain the periodicity 
of the quasiperiodic potential on the ring. Averaging 
over different /3 values is then no longer necessary since 
the computed ground-state energy of the finite-density 
Hamiltonian does not depend on f3. 

For a convenient comparison with results obtained 
for the disordered systems ||l^,|l^, we would now want 
to study the phase sensitivity at various densities, say, 
p — 1/2. However, we then have to use M = 34, 144, 610, 
. . .. For third filling, we would be restricted to M = 21, 
144, 987 and so on. In principle, such large system sizes 
M > 100 can be treated within the DMRG [Ssl, but the 
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necessary accuracy (« 10~^) is very hard to obtain for 
periodic boundary conditions. In addition, if the phase 
sensitivity decreases with system size, there is no possi- 
bihty to decide whether the decrease is algebraic or ex- 
ponential because the obtained value of MAE for such 
large system sizes is already zero within the computa- 
tional accuracy. 

For example, let us comment on p = 1/2. First, we 
check the occurrence of an MIT near /i = 1 at fixed M. 
Calculating the phase sensitivity at fixed M for different 
values of p and V, we find that there is a clear transi- 
tion near = 1, as seen in Fig. |l|. Whereas for small p 
the phase sensitivity for non-interacting particles is larger 
than for interacting — as it is in the clean case — , in- 
teraction slightly enhances the phase sensitivity in the 
extended regime for intermediate 0.5 ^ A* ^ 1. Near 
p = 1, it seems that the MIT is not shifted appreciably 
within the accuracy of our calculation. As shown in Fig. 
|l|, fitting the data results in a power-law behavior close 
to p — 1 with slightly larger pc = 1.025 for attractive 
interaction and slightly smaller pc = 0.95 for repulsive. 
Also, close to the transition, the phase sensitivity can be 
described by MAE cx {pc — pY with v k, 1. This is 
similar to the TIP situation where v = 1 was found for 
the behavior of the finite-size-scaled localization lengths 
fl@- 

In the non-interacting case, we can also study the 
length dependence of the phase sensitivity by a standard 
diagonalization routine. We investigate the system sizes 
M = 34, 144, and 610. M = 8 is excluded because 
a/(27r) = 5/8 = 0.625 differs too much from the true 
value l)/2 w 0.618. As shown in Fig. ^ the phase 

sensitivity at ^ = is different in the localized, crit- 
ical and extended regimes. Thus, for a system of free 
fermions at incommensurate density, we reproduce the 
expected transition at /ic = 1 in agreement with Refs. 
||18| . Similar plots can be made for attractive and re- 
pulsive interactions at p = 1/2 as also shown in Fig. ||. 
Unfortunately only the two system sizes M — 34 and 
144 are available due to numerical instabilities for the 
interacting system. Therefore, further conclusions ap- 
pear rather speculative. Nevertheless, we conjecture that 
the MIT at /ic ~ 1 for incommensurate densities is only 
weakly influenced by attractive and repulsive interaction 
in agreement with the results obtained in the TIP case 

IV. COMMENSURATE PARTICLE DENSITIES 

For the reasons outlined in the last section, we shall 
now turn our attention to the behavior at the commen- 
surate densities pi « hm„^oo ^n-i/^n ~ 0.618, 0.382, 
0.236, and 0.146 corresponding to i = 1, . . . , 4 in the fol- 
lowing, where a — 2npi. 



A. The Peierls-like metal-insulator-transition 

Let us first comment on the densities pi and p2 and 
note that the case of pi is identical to p2 due to the 
particle- hole symmetry in Eq. (|l]). Furthermore, the elec- 
tronic Fermi wave-vector is fcp = Trpi and therefore we 
have the resonance condition a — 2kp for the wave- vector 
a of the quasiperiodic potential. Consequently, we expect 
critical behavior as found for the periodic potential at the 
Peierls transition p^ , i.e., a transition from insulating 
to metallic phase at V = —V2 and infinite small poten- 
tial strength. Peierls Q considered about 50 years ago 
the behavior of one-dimensional tight-binding electrons 
due to a periodic lattice distortion and found the same 
resonance condition. We therefore call the transition ob- 
served in the present quasiperiodic system a Peierls-like 
transition. The critical behavior for fermions in a pe- 
riodic potential is equivalent to the critical behavior of 
the periodic hopping model. There the site-dependent 
hopping amplitude is given hy tm = 1 — ^ cos qm. In 
this so-called Su-Schrieffer-Heeger model the MIT 
occurs at V = — V2 and (5 — )■ 0, if filling factor p and 
wave vector q of the periodic hopping are commensu- 
rate, i.e., p = q/{2Tr) or q = 2kp jio). This happens for 
p — 1/2 and dimerization, one-third filling and trimer- 
ization, and so on [Q. For increasing (5, Vmit decreases. 
The quasiperiodic model is expected to show this tran- 
sition, at densities pi and p2, because a = 2kF in these 
cases. We have checked numerically that the Peierls-like 
transition is indeed observed in these cases. 

Even for the other commensurate densities ps, /04, . . ., 
where the resonance condition is not strictly fulfilled, the 
critical behavior is still expected to be similar to the 
Peierls-case as shown in Ref. In our case, if p is 

larger than a certain minimal prain value, we find that 
the phase diagrams for p^ and p4 are dominated by lo- 
calized states as shown in Fig. ^. This localized regime 
extends further on to larger interactions, —1.4 ^V< 00, 
for p > Q. However, there also exists a sizeable region of 
extended states starting roughly at F = — 1.4 for small p 
and extending until V « —2.5 for larger p. As an exam- 
ple how well the transition point is defined numerically 
consider Fig. ^ where we plot the phase sensitivity ver- 
sus interaction strength. The localization for V > —1.7 

— where MAE decreases with system size — is caused 
by the quasiperiodic potential and for V < —2 by the 
clustering of the fermions in the phase-separated state. 
In between, the quasiperiodic potential is irrelevant and 
the ground-state wave-function is extended as for p = 

— MAE increases with system size. 

In the inset of Fig. ^ we show the transition at p — 0.4 
in detail. A similar behavior close to the transition has 
also been found in data obtained recently for the TIP case 
Jl8| , p^ , albeit for varying p and constant interaction and 
not for varying V with constant p. We may nevertheless 
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perform a finite-size-scaling analysis. A simple ansatz as- 
suming a power-law divergence of the correlation length 
^ ~ |y — V'cl"'^ at the transition yields reasonable scaling 
results as shown in Fig. ^. Taking into account non- 
linearities in the finite-size-scaling ansatz as in Ref. , 
we can determine Vc and for various values of ^ as 
shown in Fig. ^. Our results indicate that v rises from 
~ 1 at /i = to values close at 2 at = 0.5 and then drops 
back to the expected = 1 at the Aubry- Andre transition 
at ^c- We emphasize that our data, due to the limited 
system sizes available, are not sufficient to distinguish 
between the power-law behavior of a second order phase 
transition and the expected Kosterlitz-Thouless behavior 
1^ at a Peierls-like transition. 

Thus, for repulsive and weak attractive interactions — 
including the non-interacting case — , we find that the 
ground state for /j, > /imin is localized. This is in agree- 
ment with previous studies for disordered and periodi- 
cally disturbed systems . We emphasize that 
an increase of the localization lengths as predicted by 
the arguments for TIP |^ is most likely too small |l3| to 
be detected by the present accuracy. For more strongly 
attractive interactions at y « — 1.4, the situation is more 
interesting. For all densities pi and fi > /imin, the sys- 
tem shows Peierls-like behavior, i.e., a transition from 
insulating to metallic phase at 1/ « —V^- Changing the 
density, the transition occurs for decreasing particle den- 
sity at increasing fi. This observation is confirmed by 
an investigation of the energy spectrum reported in Ref. 
p2| . In the non-interacting case it contains no gap at 

= but M — 1 gaps at /ic = 1 for the given a . As 
seen in Ref. most of the gaps open successively for 
increasing fj,. Especially, the first gap at k = a/2 opens 
for /i — !■ 0. Tuning the filling factor, the Fermi points 
fall into the additional gaps in the spectrum, leading to 
insulating behavior. 



B. Comparison with the renormalization group 
treatment 

In addition, the behavior seen numerically is in agree- 
ment with the weak-coupling RG treatment of spin- 
less fermions on a Fibonacci lattice. The relevant RG 
equation shows the same critical behavior as for the 
Peierls model, namely 



MAE 



-M/i 



with ^ 



2/(2~K) 



(4) 



dime ^'^/^ 



din A/ 



= (2-i^)/i, 



(3) 



where K = tt/2 arccos {—V/2) is the Luttinger parameter 
for half filling. Since the Luttinger parameter does not 
depend strongly on the filling, it is valid to use this an- 
alytical expression for other fillings, too. In accordance 
with this RG equation we expect, that the phase sensi- 
tivity — as in a disordered system — decreases in 
the localized regime such that 



where it is assumed that the localization length ^ is the 
only relevant length scale. We emphasize that this scal- 
ing is only reasonable in the localized regime. It is clearly 
found in the numerical data for, e.g., and V — —0.6 as 
shown in Fig. |^.). Near the phase transition, it is numer- 
ically difficult to distinguish between localized (exponen- 
tial decay with system size) and critical (algebraic decay 
with system size) phases. For small system sizes and 
small /i, the decrease in the phase sensitivity is always 
algebraic, according to the RG equation. The exponen- 
tial decrease with f sets in far from the transition points 
as shown in Fig. ^). There, the exponential scaling is 
found for p > 0.5. We note, that despite the straight line 
was obtained by a fit only to the data for p — 0.7 and 
p = 0.8 the data for p = 0.5, and p — 0.6 fall likewise 
on this line. To show the deviations from this straight 
line for small p clearly, which is the most important fea- 
ture in this plot, we do not show the data for p = 0.8. 
We call this intermediate region, where no universal alge- 
braic or exponential decrease is found, transition region. 
A finite-size scaling corresponding to the previous section 
near the transition with ^ ~ \p — pd'^ does not lead to 
conclusive results. 



C. Enhancement near the first-order transition 

As already mentioned in [ p0| , the phase sensitivity 
shows a large enhancement around V = —2. We find 
that the first order phase transition at ^ = — 2 remains 
unaffected for small p ^ 0.2. However, for larger p val- 
ues, the metallic phase extends towards more negative V 
values. For p > 0.8 at Vp ^ —2.2 the phase sensitivity 
shows an unexpected sharp maximum in the delocalized 
regime, as shown in Fig. ^. A similar maximum is also 
seen for p > 1 and small system sizes (M — 34, 55), 
but the system is insulating. The exponential scaling in 
the localized region as shown in Fig. does not change 
for V > —2 and V < Vp < —2 in the case of large p. 
Therefore, we conclude that this regime belongs to the 
localized regime of the quasiperiodic potential. On the 
other hand, no scaling behavior is found for V < —2 
and small and intermediate p, lower left corner in the 
phase diagrams, Fig. ^. Here, the localization is due to 
the formation of particle clusters in the phase-separated 
state. In the crossover regime between these two local- 
ized phases a metallic state is recovered. The maximum 
in the insulating regime is a rest of this crossover. At 
intermediate and strong p and V < Vp < —2 the phase 
sensitivity drops very rapidly to zero within the numer- 
ical accuracy. This behavior is comparable to the be- 
havior at V — —2 and small p as seen in Fig. Q at the 
first order transition. Thus, it seems that the first order 
transition is moved to smaller V with increasing p. The 
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subtle interplay between the clustering and the quasiperi- 
odic potential seems to delocalize the ground-state wave- 
function in a small parameter region. However, it is nu- 
merically difficult to follow the border between these two 
completely different localized phases for stronger fi and 
— 2 < y < — 3 in our approach. 

D. Phase diagram 

In Fig. ^, we summarize our results by showing two 
phase diagrams of model for varying quasiperiodic 
potential strength /i and nearest-neighbor interaction V. 
They were obtained by studying the system size behavior 
of MAE up to M = 144 for pa = 0.236 and p4 = 0.146. 
The localized parts of the phase diagrams for weak at- 
tractive and repulsive interactions are not shown, be- 
cause they have no additional structure. For p^, the 
intermediate phase is difficult to detect and is not sepa- 
rately marked in Fig. |^. There we can only distinguish 
whether the phase sensitivity decreases or not. As clearly 
seen, pmin = for and /imin = 0.1 for p4. In sum- 
mary, the delocalized phase becomes larger for weak fi 
and smaller for strong p. with decreasing density. First, 
/^min(/04) > P'minips) and sccoud, the crossover region 
beyond V = —2 becomes smaller and extends only to 
p = 0.9. 

V. CONCLUSIONS 



potential strength. Thus, although we find a strong en- 
hancement of the phase sensitivity at the first order tran- 
sition, the physics of the model at commensurate den- 
sities is dominated by the Peierls resonance condition 
which becomes irrelevant only for strong attractive in- 
teraction. A remainder of the single-particle localization 
is still found in this case — the region of the phase dia- 
gram, where the ground-state wave-function is extended, 
extends at most to /i = 1. In addition, we have found 
a change in the exponent of the localization length from 
i/w2toj'«lat^=l. In summary, we have seen 
that the single-particle localization in the Aubry- Andre 
model is not influenced by interaction. At commensurate 
densities, other effects such as a Peierls-like commensu- 
rability become important and dominate the transport 
properties. 
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In this work, we have studied the influence of interac- 
tions on an metal insulator transition in one dimension 
using the DMRG method. Thereby, we have compared 
the interplay of disorder and interactions for the quantum 
system (|l|) in a quasiperiodic potential at incommensu- 
rate and commensurate densities. Our results suggest 
that the dclocalization found previously for low density 
TIP in the localized phase cannot simply be extrapo- 
lated to the finite-density situation. In the case of in- 
commensurate densities with nearest-neighbor (repulsive 
or attractive) interaction the metal insulator transition 
is found at /ic ~ 1. The critical exponent, v ^ I, is 
in agreement with the results of transfer-matrix-method 
calculations and finite-size scaling This indicates 

that the influence of interaction at this particular type of 
metal insulator transition is not strong enough to change 
the universality class of the Aubry- Andre model. 

For commensurate densities p^ and p4, we have 
deduced the phase diagrams of system for varying 
quasiperiodic potential strength p and interaction V. 
The numerically accessible filling factors are commensu- 
rate with the quasiperiodic potential and yield a Peierls- 
like behavior of the system with a metal insulator tran- 
sition at attractive interaction and small quasiperiodic 
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FIG. 1. Phase sensitivity versus ^ at fixed system size 
M — 34 for V = —1.5 (squares), (circles), 1.5 (diamonds) at 
(incommensurate) half filling. The lines are fits corresponding 
MAE ~ l/i - pc|. 
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FIG. 2. Phase sensitivity versus system size for V — ~1 
(squares), (circles), 1 (diamonds) at (incommensurate) half 
filling. The three cases n = 0.9 (solid lines) , fj, = 1 (dashed 
lines), and fj, = 1.1 (dotted lines) are compared. 
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FIG. 3. Phase diagrams of the system described by Eq. 
in terms of quasiperiodic potential strength ^ and interactioi 
V for p = p3 (left side) and p — p4 (right side). An extendec 
ground state wave function is marked by •, a locahzed by o 
In addition, the transition regime (see text) is marked witi 
a shaded □ in the left figure. The solid lines indicate th( 
first-order transition at V = — 2 and p = 0. 
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FIG. 4. Phase sensitivity versus interaction for three dif- 
ferent system sizes, density p3 and p — 0.4. The lines for 
M = 34 and 89 are guides to the eye, highlighting the differ- 
ent finite-size-scaling behaviors in the various regimes. The 
dotted lines indicate the transition between the localized and 
the extended wave functions. Inset: Additional data close to 
the transition and M = 144 (A) included. The lines denote 
the curves constructed by non-linear finite-size scaling. 




FIG. 5. Scaling function (solid line) and scaled data points 
in the interval V € [—1.8, —1.4] close to the transition for 
p — 0.4 and system sizes M — 34,55,89, 144. Inset: Scaling 
parameter ^ as a function of V. 
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FIG. 6. Estimates of K (o) and v (□) obtained for 
P3 by non-linear finite-size scaling at various p. The open 
symbols denote various fit functions and initial parameters, 
whereas the closed symbols indicate the best (according to 
X^-statistics) of all these fit for a given value of p. Error bars 
mark the errors resulting from the Levenberg-Marquardt fit- 
ting method. 
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FIG. 7. (a) Phase sensitivity versus scaled system size 
^j^2/(2-if)^ where M = 34, 55, and 89, for V = -0.6, filling 
p = p3, and various potential strengths. The line indicates a 
plot of Eq. (§) with K = 7r/2 arccos(0.3) « 1.24. (b) Phase 
sensitivity versus scaled system size, M — 34, 55, and 89, for 
V = —1.5 and p^. The straight line is a fit to the data for 
p = 0.7 and 0.8. 
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FIG. 8. Phase sensitivity versus interaction at ps for dif- 
ferent fj, at size M — 34 (small grey symbols) and AI = 55 
(large open symbols). We compare two system sizes to indi- 
cate the regimes with the localized and extended ground state 
wave function. The lines connect the data points for different 
fj, and Af values. 
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